AL609874 


VUF-1013 

FINAL  REPORT 


VELA  UNIFORM 


$HML 


SPONSORED  BY  THE  ADVANCED  RESEARCH  PROJECTS  AGENCY  OF  THE 
DEPARTMENT  OF  DEFENSE  AND  THE  U.S.  ATOMIC  ENERGY  COMMISSION 

COPY  — T^OF - '  '  ' 

NEVADA 


1963 


ANALYSIS  OF  SHOAL  DATA  ON  GROUND 
MOTION  AND  CONTAINMENT 

Roland  F.  B^ers,  Inc. 

December  4,  1964 


Issuance  Date:  January  11,  19(15 


LEGAL  NOTICE 

This  report  was  prepared  as  an  account  of  Government  sponsored  work.  Neither  the  United 
SL'iles.  nor  the  Commi*«i'>n.  »^nr  any  pr»rsor.  acting  on  behaU  ''f  the  Commission; 

A.  Makes  any  warranty  or  represenUlion.  expressed  or  Implied,  with  respect  to  the  accu¬ 
racy.  completeness,  or  usefulness  of  the  information  contained  in  this  report,  or  that  the  use 
of  any  Information,  apparatus,  method,  or  process  disclosed  In  this  report  may  not  infringe 
privately  owTied  rights:  or 

B.  Assumes  any  liabilities  with  respect  to  the  use  of,  or  for  damages  resulting  from  the 
use  of  any  Information,  apparatus,  method,  or  process  disclosed  in  this  report. 

As  used  In  the  above,  “person  acting  on  behalf  of  the  Commission'*  includes  any  em¬ 
ployee  or  contractor  of  the  Commission,  or  employee  of  such  contractor,  to  the  extent  that 
such  employee  or  contractor  of  the  Commission,  or  employee  of  such  contractor  prepares, 
disseminates,  or  provides  access  to,  any  Information  pursuant  to  his  employment  or  contract 
with  the  Commission,  or  hla  employment  with  such  contractor. 


This  report  has  been  reproduced  directly  from  the  best 
available  copy. 

Printed  in  USA.  Price  $4.00.  Available  from  the  Clearing¬ 
house  for  Federal  Scientific  and  Technical  Information,  Na¬ 
tional  Bureau  of  Standards,  U.  S.  Department  of  Commerce, 
Springfield,  Va. 


Reproduced  From 

Best  Available  copy 


NUCLEAR  EXPLOSIONS — PEACEFUL  APPLICATIONS 

TID-4S0O(:nth  Ed.) 


PR0.7ECT  SHOAL 


Final  Report 


ANALYSIS  OF  SHOAL  DATA 
ON  GROUND  MOTION  AND  CONTAINMENT 

SAFETY  PROGRAM 


December  4,  1964 


Roland  F.  Beers,  Inc. 
Alexandria,  Virginia 


CONTENTS 

Page 

ABSTRACT  1 

CHAPTER  I  INTRODUCTION  3 

1.1  Background  3 

1.2  Objectives  6 

1.3  Ir.struraentation  7 

1.4  Theory  9 

CHAPTER  II  PREDICTIONS  13 

2.1  Geologic  Environment  13 

2.2  Phenomenology  and  Containment  15 

2.3  Seismic  Effects  35 

CHAPTER  III  ANALYSIS  AND  INTERPRETATION  40 

3.1  Phenomenology  and  Containment  40 

3.2  Seismic  Effects  42 

CHAPTER  IV  CONCLUSIONS  AND  RECOMMENT5ATIONS  52 

APPENDIX  A  54 

APPEiTDIX  61 

TABLES 

2.1  Comparison  of  Physical  Properties  of 

Shoal  and  Climax  Stock  Granites  64 

2.2  Predictions  of  Amplitudes  of  Ground 

Motion  65 

3.1  Peak  Surface  Particle  Acceleration  in  g's  66 

3.2  Peak  Surface  Particle  Displacement  in 

Centimeters  68 

3.3  Peak  Surface  Particle  Velocity  in  cm/sec  70 

3.4  Stations  Used  in  Comparison  of  Data 

Measured  Parallel  and  Perpendicular  to 
Direction  of  Major  Faulting  72 

3.5  Ratio  of  73 

ad 

FIGURES 

1.1  Locations  of  Instrument  Stations  74 

1.2  Cross  Section  A-A'-A"  75 

1.3  Geologic  Map  76 

in 


Preceding  Page  Blank 


Cas’TENTS  (continued) 


Page 

FIGURES 

2.1  Predicted  and  Observed  Permanent  Radial 

Displacement  and  Predicted  Strain  for 

the  Hardhat  Event  77 

2.2  Predicted  Strain  and  Permanent  Radial 

Displacement  for  the  Shoal  Event  78 

2.3  Spall  Depth  as  a  Function  of  Tensile 

Strength  for  Shoal  Grajtiite  79 

2.4  Predicted  Peak  Surface  Particle  Accel¬ 

eration  as  a  Function  of  Horizontal 
Distance,  Shoal  Event  80 

2.5  Predicted  Pe2dc  Surface  Particle  Dis¬ 

placement  as  a  Function  of  Horizontal 
Distance,  Shoal  Event  81 

2.6  Predicted  Peak  Surface  Particle  Velocity 

as  a  Function  of  Horizontal  Distance, 

Shoal  Event  82 

2.7  Area  Map  Showing  Predicted  Distances  to 

Significant  Ground  Motions  83 

3.1  Vertical  Component  of  Acceleration  84 

3.2  Radial  Component  of  Acceleration  85 

3.3  Transverse  Component  of  Acceleration  86 

3.4  Resultant  Vector  of  Acceleration  87 

3.5  Vertical  Component  of  Displacement  88 

3.6  Radial  Component  of  Displacement  89 

3.7  Transverse  Component  of  Displacement  90 

3.8  Resultant  Vector  of  Displacement  91 

3.9  Vertical  Component  of  Velocity  92 

3.10  Radial  Component  of  Velocity  93 

3.11  Transverse  Component  of  Velocity  94 

3.12  Resultant  Vector  of  Velocity  95 

3.13  Comparison  of  Resultant  Vector  of  Accel¬ 

eration  Measured  Parallel  and  Perpen¬ 
dicular  to  Major  Fault  Trend  96 

3.14  Comparison  of  Resultant  Vector  of  Dis¬ 

placement  Measured  Parallel  and  Perpen¬ 
dicular  to  Major  Fault  Trend  97 

3.15  Comparison  of  Resultant  Vector  of  Veloci¬ 

ty  Measured  Parallel  and  Perpendicular 
to  Major  Fault  Trend  98 


IV 


CONTENTS  (continued) 


FIGURES 

3.16  Comparison  of  Predicted  and  Observed 

Resultant  Vector  of  Acceleration 

3.17  Comparison  of  Predicted  and  Observed 

Resultant  Vector  of  Displacement 

3.18  Comparison  of  Predicted  and  Observed 

Resultant  Vector  of  Particle  Velocity 


Page 


99 

100 

101 


REFERENCES 


102 


ABSTOACT 


Predictions  and  evaluations  of  ground  motion  and  con¬ 
tainment  were  made  for  the  12.5  kt  Shoal  event.  A  compre¬ 
hensive  ground  motion  instrumentation  program  was  dei...gned. 
This  program  included  four  linear  instrument  arrays  radi¬ 
ating  northeast,  northwest,  southeast  and  southwest  with  a 
total  of  21  stations.  Eight  additional  stations  were  em¬ 
placed  at  points  of  interest.  Of  the  29  stations,  21  were 
on  hardrock  and  eight  were  on  alluvium. 

Analysis  of  the  seismic  records  included  corrections 
for  instrument  response  and  derivation  of  non-recorded  mo¬ 
tions  by  differentiation  or  integration.  It  was  determined 
that  p>eak  motions  on  alluvium  were  greater  than  on  hardrock 
by  a  factor  of  2.1  for  acceleration,  3.7  for  displacement 
and  3,4  for  velocity.  The  rates  of  attenuation  indicated 
that  p>eak  motions  on  alluvium  were  not  due  to  surface 
waves.  The  prediction  equations  for  acceleration  and  dis¬ 
placement  were  found  to  represent  reasonable  upper  limits 
of  motion. 

The  objective  of  the  containment  evaluation  was  to  as¬ 
sure  that  appropriate  measures  had  been  taken  to  preclude 
unacceptable  venting  of  radioactive  isotope  to  the 


atmosphere.  Consideration  was  given  to  possible  venting 
both  at  the  time  of  the  explosion  and  later  when  the  cavity 
collapsed  and  a  chimney  formed.  The  adequacy  of  stemming 
of  all  man-made  openings  as  well  as  the  effect  of  geologic 
hazards  such  as  faults  or  fissures  were  evaluated  in  the 
containment  analysis.  The  possibility  of  a  secondary  ex¬ 
plosion  due  to  water  entering  the  cavity  was  also  examined. 
It  was  concluded  that  Shoal  would  be  contained  to  the  same 
degree  as  Kardhat  and  that  the  stemming  plan  was  adequate. 

The  prediction  of  containment  was  substantiated  by 
the  test  results.  Limited  post-shot  exploration  allows 
the  following  comparisons:  cavity  radius  -  24.7  meters 
predicted,  25.6  0.6  meters  observed;  chimney  height  - 

110  meters  predicted,  108  meters  observed;  depth  of  spall¬ 
ing  -  122  metex^s  predicted,  120+  meters  indicated;  radius 
of  spalling  -  <914  meters  predicted,  794  meters  observed. 
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CHAPTER  I 


INTRODUCTION 

1 . 1  BACKGROUND 

Rcland  F.  Beers,  Inc,,  under  Contract  AT(  29-2 )- 1163 
with  the  Nevada  Operations  Office  (NVOO) ,  U.  S.  Atomic 
Energy  Commission  (AEC)  is  responsible  to  the  Operational 
Safety  Division  of  that  office  for  the  following: 

a)  Prediction  of  ground  motions 

b)  Evaluation  of  containment  of  underground  nuclear 
explos ions 

The  12.5  kt  Shoal  event  was  the  second,  and  largest, 
contained  underground  nuclear  test  conducted  in  granite. 
This  event  presented  special  prediction  problems  due  to 
the  absence  of  directly  applicable  experience  and  its  lo¬ 
cation  off  the  Nevada  Test  Site  in  the  Sand  Springs  moun¬ 
tain  range  near  Fallon,  Nevada. 

Ground  motion  resulting  from  underground  testing  of 
nuclear  devices  is  a  substantial  as  well  as  continuing 
operational  safetv  problem.  Current  methods  of  predict¬ 
ing  ground  motion  employ  empirical  relationships  derived 
from  observed  data.  Proper  use  of  these  relationships 
for  prediction  purposes  requires  that  the  test  conditions 
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for  which  predictions  are  to  be  made  be  coinp)arable  with 
those  which  provided  the  data  used  in  deriving  the  rela¬ 
tionships  . 

Scaling  of  ground  motion  for  yield  and  distance  are 
fundamental  steps  in  making  predictions.  To  perform  such 
scaling,  the  variation  in  amplitude  with  yield  and  dis¬ 
tance  must  be  defined.  Appreciable  departure  in  yield, 
source  geometry,  method  of  emplacement  and  physical  prop¬ 
erties  of  the  transmission  path  from  previous  experiments 
substantially  affect  manner  of  scaling  and  the  agreement 
between  predicted  and  observed  ground  motions. 

The  results  of  the  earlier  5.9  kt  Hardhat  event  in 
granite  on  NTS  were  not  used.  The  limited  span  of  dis¬ 
tances  at  which  ground  motion  data  were  recorded  from 
Hardhat  and  the  high  degree  of  scatter  in  sur f ace -recorded 
data  resulted  in  an  ambiguous  rate  of  attenuation  of  ground 
motion  with  distance.  Additionally,  as  Hardhat  was  the 
only  test  conducted  in  granite  prior  to  Shoa] ,  the  empiri¬ 
cal  scaling  of  amplitude  with  yield  was  indeterminate. 

To  find  an  alternate  basis  for  making  predictions,  the  dis¬ 
placement  potentials  were  compared  for  the  different  media 
in  which  tests  have  been  conducted  to  date.  This 
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comparison  suggested  that  potentially  more  accurate  re¬ 
sults,  especially  at  larger  distances,  should  be  obtained 
using  tuff  as  a  basis  for  predicting  ground  mentions.  Ac¬ 
cordingly,  prediction  relationships  were  deriveci  from  the 
results  of  the  Hardtack-Plumbob  experiments  in  Wainier 
Mesa  tuff.  However,  differences  between  the  geology  of  the 
Shoal  site  and  Rainier  Mesa  introduced  substantial  uncer¬ 
tainty.  The  scope  of  the  transmission  problem  is  defined 
in  Section  1.4. 

The  possibility  of  a  containment  failure  is  a  re¬ 
curring  problem.  This  hazard  becomes  particularly  impor¬ 
tant  when  an  off-site  event  is  considered.  In  analyzing 
the  Shoal  event  from  this  viewpoint,  the  Hardhat  event 
provided  the  only  previous  exp>erience  in  a  granitic  med  .m 
upon  which  to  base  predictions.  The  prediction  of  pheno¬ 
menology,  therefore,  relied  heavily  upon  the  available 
Hardhat  data.  Where  conditions  were  dissimilar,  theoreti¬ 
cal  analysis  and  extrapolation  of  experience  from  other 
media  were  employed.  These  techniques  permitted  an  evalua¬ 
tion  to  be  made  of  unique  conditions  such  as  the  stemming 
design  and  the  presence  of  a  significant  fault  in  the  vi¬ 
cinity  of  the  weapon  point. 
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1.2  OBJECTIVES 


In  regard  to  contalnoent ,  the  purpose  was  to  provide 
the  Operational  Safety  Division  with  an  independent  evalua¬ 
tion  of  the  probability  of  containment  and  the  adequacy  of 
the  stemming  concept  and  plan.  In  {performing  this  function, 
dependence  was  placed  on  theoretical  considerations  and  on 
empirical  data  collected  in  large  measure  by  others.  This 
supply  of  data  was  arranged  primarily  by  the  NVOO,  but  al¬ 
so  wa*5  obtained  by  direct  contact  with  the  Technical  Direc¬ 
tor  where  appropriate. 

To  evaluate  the  effect  of  the  varied  geology  on  ground 
motion,  a  comprehensive  ground  motion  measurement  program 
was  designed.  Instrumentation  program  design  objectives 
for  the  Shoal  event  were  to  obtain  ground  motion  data  per¬ 
mitting  evaluation  of  the  following; 

1)  The  variation  in  the  amplitude  of  ground  motions 
on  granite  and  alluvium  with  distance. 

2)  The  effect  of  faulting  on  ground  motion. 

3)  The  effect  of  thickness  of  alluvium  on  the  am¬ 
plitude  and  attenuation  of  amplitude  with  dis¬ 
tance  . 

4)  The  partition  of  compress ional  wave  energy  into 
surface  waves  at  a  lithologic  discontinuity. 
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5)  The  directional  response  characteristics  of 
the  shot  medium. 

6)  The  degree  of  scatter  in  recorded  ground  mo¬ 
tions  . 

7)  The  ratio  of  (u^/ad)  where  ”u"  is  particle  ve¬ 
locity,  "a"  is  p>article  acceleration,  and  "d” 
is  particle  displacement. 

1.3  INSTRUMENTATION 

To  accomplish  the  objectives  listed  in  Section  1.2, 

21  recording  stations  were  deployed  in  four  linear  arrays 
radiating  northeast,  northwest,  southeast  and  southwest 
from  ground  zero  (see  Figure  1.1).  Each  of  these  21  sta¬ 
tions  had  aui  instrument  complement  of  three  component  ac¬ 
celerometers  aoid  three  component  displacement  meters. 

Three  of  the  21  recording  stations  also  included  instruments 
for  recording  three  components  of  velocity.  These  instru¬ 
ments  were  located  on  the  southeast  line.  Besides  the  above, 
vertical  component  velocity  measurements  were  made  at  five 
additional  stations  at  disteuice.  Each  of  the  vertical  com¬ 
ponent  measurements  represents  the  output  of  a  circular  ar¬ 
ray  of  six  vertical  velocity  instruments  emplaced  at  widely 
separated  points  surrounding  the  test  site.  An  additional 
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three  vibration  meters  measuring  radial  and  transverse  dis¬ 
placement  were  located  at  the  Nevada  Scheelite  Mine,  Gabbs 
Mine  and  the  Naval  Air  Station. 

Of  the  21  recording  stations  deployed  in  arrays,  four¬ 
teen  were  located  on  gramite  and  seven  were  located  on  al¬ 
luvium.  To  evaluate  the  transmission  properties  of  gran¬ 
ite,  the  positions  of  the  recording  arrays  located  on  gran¬ 
ite  were  given  preferential  azimuthal  orientations  related 
to  the  more  predominant  NE-SW  system  of  faulting.  The  ar¬ 
rays  oriented  northeast-southwest  essentially  parallel  this 
system  of  faulting  while  the  arrays  oriented  northwest- 
southeast  intersect  this  faulting  essentially  at  right  an¬ 
gles.  Evaluation  of  the  effect  of  different  lithology  was 
provided  for  by  extending  the  arrays  from  granite  to  allu¬ 
vium.  Evaluation  of  the  effect  of  different  thicknesses  of 
alluvium  was  provided  for  by  locating  one  of  the  arrays  on 
the  thick  alluvium  of  Fairview  Valley  and  tl>e  other  on  the 
thinner  alluvium  of  Fourmile  and  Eightmile  Flats.  The  cir¬ 
cular  arrays  of  vertical  component  velocity  meters  emplaced 
at  distamce  were  on  hardrock  of  unknown  lithology.  These 
were  not  part  of  the  safety  studies  program.  The  vibration 
meters  were  emplaced  to  monitor  ground  motions  at  sites  of 
engineering  structures. 
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1.4  THEORY 


1.4.1  Seismic  Effects.  At  the  Shoal  site  consider¬ 
able  variation  in  the  amplitudes  of  the  particle  motions 
was  expected  to  occur  both  horizontally  and  directionally 
due  to  asymmetries  in  the  geometry  and  geology  of  the 
transmission  model.  Figure  1.2  is  a  cross  section  showing 
the  attitudes  of  the  granite-alluvium  contact  on  either 
side  of  the  site.  The  eastern  contact  of  the  granite  with 
the  alluvium  is  abrupt  with  nearly  vertical  dip  and  with  a 
great  depth  of  alluvium  directlv  adjacent  to  the  granite. 
The  western  edge  of  the  granite  shows  much  less  dip  with 
alluvium  forming  a  wedge  over  the  granite. 

The  geologic  map  of  Figure  1.3  shows  the  trend  of 
faulting  at  the  Shoal  site  as  far  as  it  is  known.  The 
major  trend  is  northeast-southwest  on  the  basis  of  surface 
evidence  but  it  should  be  noted  that  a  major  fault  trend¬ 
ing  east-west  lies  at  a  radial  distance  of  approximately 
34  meters  north  of  the  working  p>oint.  The  multiplicity 
of  southeast-northwest  trending  features  indicated  on  the 
map  are  dikes. 

The  Shoal  model  described  above  departs  from  the  homo¬ 
geneous,  isotropic  conditions  usually  assumed.  The  many 
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^faults,  Joints  and  dikes  in  the  granite  were  (xpected  tf) 
have  an  effect  on  azimuthal  variation  of  ground  motion. 

The  energy  transmission  parallel  to  the  northeast-southwest 
set  would  be  enhanced  as  compared  to  energy  transmission  at 
right  angles  to  the  fault  set.  This  means  that  energy 
reaching  the  alluvium  of  Fourmile  Flat  and  Fairview  Valley 
at  the  points  closest  to  the  working  point  may  be  less 
than  that  normally  expected.  The  energy  along  the  Sand 
Springs  Range  would  be  about  as  expected  to  the  south  and 
slightly  lower  than  expected  to  the  north  due  to  the  large 
cast-west  fault  Just  north  of  the  working  point . 

The  second  aspect  of  the  model  relates  to  the  pres¬ 
ence  of  the  alluvium  in  the  adjacent  Fourmile  Flat  and 
Fairview  Valley  areas.  The  question  here  was  whether  sur¬ 
face  waves  generated  at  the  contact  between  granite  and  al¬ 
luvium  would  produce  larger  ground  motions  at  distances  on 
the  alluvium  than  those  resulting  from  compress i  onai  (body) 
wave  arrivals.  The  amplitude  of  the  ground  motion  due  to 
surface  waves  for  a  two-media  model  depends  on  diverse  fac¬ 
tors.  Some  of  these  are:  the  dista-nce  of  the  granite- 
alluvium  contact  from  the  source,  the  dip  of  the  contact, 
the  partition  of  energy  from  corapressional  waves  into  sur¬ 
face  waves,  and  the  thickness  of  alluvium  at  the  contact 
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and  along  the  transmission  path  to  the  point  of  interest. 
Although  Ewing,  Jardetzky  and  Press  (Reference  19)  have 
developed  surface  wave  theory  for  specialized  cases,  there 
was  no  known  theory  applicable  to  the  Shoal  model.  Together 
with  this,  there  was  also  a  lack  of  empirical  data  on  par¬ 
tition  of  energy  from  compressional  into  surface  waves. 

On  the  other  hand,  to  predict  the  ground  motions  at 
the  surface  of  the  alluvium  due  to  compressional  waves,  it 
is  required  to  increase  the  amplitudes  of  predicted  free- 
field  particle  motions  by  some  factor.  This  factor  proba¬ 
bly  depends  on  the  interrelationship  of  wave  p^eriod  to 
bed  thickness.  The  factor  two  has  been  commonly  used  in 
the  past  to  convert  from  surface  to  subsurface  motion, 
but  this  does  not  include  effects  of  transmission  from  one 
rock  type  to  another.  Because  of  these  uncertainties,  no 
factor  was  applied  in  making  pnredictions  of  Shoal  seismic 
effects  on  alluvium. 

1.4.2  Phenomenology  and  Containment.  The  evaluation 
of  containment  of  the  Shoal  event  was  based  on  procedures 
developed  as  part  of  a  long  range  program  directed  to  the 
continuing  improvement  of  confidence  in  predictions. 
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Currently,  greater  dependence  is  placed  on  empirical  rela¬ 
tionships,  and  it  is  probable  that  experience  will  always 
be  a  major  factor  in  evaluation.  However,  a  large  part  of 
the  improvement  will  come  from  greater  knowledge  of  the 
pressure- time  history  aund  temper ature- 1 ime  history  of  un¬ 
derground  explosions.  This  is  being  actively  pursued, 
and  some  of  the  results  were  employed  in  the  development 
of  predictions  and  the  evaluation  for  Shoal. 
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CHAPTER  II 


PREDICTIONS 

2.1  GEOLOGIC  ENVIRONMENT 

The  local  geologic  and  topographic  setting  of  the 
Sand  Springs  Range  and  the  surrounding  terrain  has  been 
amply  provided  by  the  text,  maps  and  cross  sections  of 
Reference  1.  A  cross  section  and  a  portion  of  a  map 
from  this  reference  are  included  in  Chapter  1. 

The  lithology  of  the  granitic  shot  medium  is  also 
described  in  Reference  1.  Its  physical  properties  have 
been  compiled  from  several  sources  and  are  tabulated  in 
Table  2.1. 

The  significant  structural  features  of  the  region 
include  a  well-developed  fracture  cleavage  which  per¬ 
sists  throughout  the  granitic  medium.  Generally,  planes 
are  ,25  to  .51  cm  apart,  trending  northeast,  with  a  near 
vertical  dip.  Also,  steeply  dipping  wel 1 -developed 
joints  trend  northwesterly  (Reference  1). 

The  entire  Sand  Springs  Range  is  heavily  faulted, 
the  strongest  trend  being  northeasterly.  However,  geo¬ 
logic  evidence  from  ECH-D  (Reference  1),  the  cast  drift 
(Reference  2),  and  test  holes  B-1  through  B-4,  (Refer¬ 
ences  3  and  4)  indicate  with  certainty  that  no  signifi¬ 
cant  faults  exist  within  a  30  meter  radius  of  the 
Working  Point  (W.P,).  The  test  holes  indicate  that  the 
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rtearest  fault  is  the  so-called  Fault  E  which,  prior  to 
detonation,  lay  at  a  radial  distance  of  about  34  meters 
from  the  W.P,  Correlation  with  its  probable  surface 
exposure  indicates  that  the  fault  strikes  approximately 
east-west  and  dips  70  degrees  south.  Its  true  pre-shot 
thickness,  where  intersected  by  ECH-D,  was  about  34 
meters , 

Because  of  considerable  core  loss  in  drilj^ing  ECH-D 
in  the  fault  zone,  the  material  character  was  not  directly 
observed.  However,  geophysical  logs  of  this  hole  suggest 
the  zone  is  composed  of  slightly-to-completely  decomposed 
granite  containing  water,  most  of  which  is  fixed  mole- 
cularly  by  hydration  of  feldspar,.  Large  voids  are  prob¬ 
ably  absent  and  the  fault  material  is  considered  to  be 
mostly  impervious. 

The  next  nearest  known  fault  is  Zone  C,  originally 
identified  on  the  surface  and  later  intersected  by  the 
east  drift  at  Static  n  8  42 .  It  strikes  approximately 

north  -  22  degrees  -  east  and  dips  about  74  degrees  south¬ 
east.  Its  true  thickness  is  about  3.65  meters,  and  is 
composed  chiefly  of  slickensided  gouge  (Reference  2). 

Ground  water  within  the  granite  occurs  within  joints, 
fractures  and  faults  which  are  probably  incompletely  con¬ 
nected  and  which  may  or  may  not  share  a  vague  and  irregular 
water  table.  Initially,  such  a  surface  was  found  at  a 
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depth  of  about  294  meters  in  Hole  ECH-D  (Reference  1). 

A  later  pump  test  prior  to  underground  excavation  in¬ 
dicated  the  water  level  may  have  dropped  to  a  depth  of 
316  meters.  It  seems  likely  that  many  of  the  local 
water-bearing  zones  have  been  lowered  appreciably  as  a 
result  of  drainage  by  subsequent  underground  excavation. 
Nonetheless,  the  continued  heavy  flow  from  Hole  E-1 
(Reference  5)  indicated  that  considerable  ground  water 
with  appreciable  head  remained  in  the  vicinity  at  the 
shot  time.  Furthermore,  the  discontinuous  nature  of  the 
ground  water  indicates  that  a  water-bearing  zone  “;ould 
remain  essentially  unchanged  by  drift  excavation  and 
could  retain  its  initial  hydrostatic  head.  Therefore, 
as  stated  in  Reference  6,  it  is  conceivable  that  such 
a  water-bearing  zone  exists  within  the  fracture  limit 
of  the  detonation. 

C 

2.2  PHENOMENOLOGY  AND  CONTAINMENT 

The  following  sections  provide  the  basis  for  the 
predictions  shown  in  the  abstract  pertaining  to  close- 
in  phenomena  and  containment  problems. 

2.2,1  Radius  of  Cavity.  The  predicted  cavity 
radius  was  24.7  meters.  this  was  based  on  cube-root 
scaling  from  Hardhat  where  the  observed  cavity  radius 
was  19.2  meters  (Reference  7).  The  yield  used  for 
Shoal  was  the  measured  12.5  kt  and  that  for  Hardhat  was 
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5.9  kt.  The  chamber  sizes  for  Hardhat  and  Shoal  are 
comparable.  Accordingly,  cube-root  scaling  should  be 
applicable.  A  check  was  performed  establishing  the 
peak  pressure  at  the  cavity  radius  for  F^ardhat  and  coa- 
puting  the  distance  for  equal  pressure  for  Shoal.  This 
confirmed  the  cube -root  scaling  procedure. 

2,2.2  Cracking  Radius.  I?adial  cxacking  for  the 
Shoal  event  was  predicted  to  extend  to  approximately 
159  meters  in  the  free  field  and  possibly  greater  dis¬ 
tances  where  the  presence  of  openings  produces  a  free- 
face  effect . 

The  extent  to  which  fractures  will  be  propagated 
due  to  a  subsurface  detonation  has  been  found  to  depend 
upon  the  physical  properties  of  the  shot  medium.  It  is 
therefore  necessary  to  rely  upon  past  experience  in  the 
medium  and  any  experimental  data  which  have  a  bearing  on 
the  problem.  From  Hardhat  it  was  found  that  post-shot 
fractures,  observable  in  the  tunnel,  extended  to  a  dis¬ 
tance  of  159  meters  from  the  W.P.  (Reference  8);  free- 
field  fractures  would  be  expected  to  be  less  but  no 
maximum  value  is  recorded  for  these.  Radioactivity  was 
encountered  in  the  end  of  the  Hardhat  drift  and  there¬ 
fore  continuous  channels  for  transport  must  have  existed 
out  to  this  distance,  a  radius  of  approximately  61  meters. 
This  may  then  constitute  a  minlnuin  free-field  fracture 
radius . 
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The  above  values  were  used  to  predict  tlie  extent 
of  radial  fracturing  anticipated  for  the  Shoal  event  by 
considering  the  permanent  radial  displacements  and  as¬ 
sociated  strains  for  the  two  events.  The  predicted  and 
observed  radial  displacements  were  derived  by  a  method 
similar  to  that  outlined  in  Appendix  A  and  are  shown  in 
Figures  2.1  and  2.2.  The  use  of  strain  data  in  obtaining 
predictions  of  fracturing  radius  is  logical  if  the  mode 
of  deformation  is  considered.  As  the  cavity  expands,  the 
forces  existing  in  the  free  field  are  a  radial  compression 
coupled  with  a  tangential  expansion.  This  situation  may 
be  approximated  by  uniaxial  compression  tests  since  the 
stresses  created  by  tangential  tension  cannot  become  large 
for  geologic  materials.  The  permanent  strain  of  159  meters 
for  Hardhat  was  0.125  percent  as  obtained  from  Figure  2.1. 
Using  this  strain,  fracturing  would  extend  to  a  distance 
of  201  meters  for  Shoal.  Similarly,  a  minimum  value  was 
established  for  Shoal  by  considering  the  strain  at  61  meters 
for  Hardhat  since  radioactive  contamination  entered  the 
drift  at  this  radius.  This  yielded  a  radial  distance  of 
about  76  meters  for  an  equivalent  strain  at  Shoal . 

Since  the  maximum  fracturing  radius  observed  in  tun¬ 
nels  is  certainly  greater  than  in  the  free-fleld  condition, 
it  is  necessary  to  obtain  estimates  of  the  latter. 

From  results  of  three  uniaxial  compression  tests  on 
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the  Shoal  granite  cores  by  the  Waterways  Experiment  Sta¬ 
tion,  an  average  uaximuio  compressive  strain  of  0.25  per¬ 
cent  was  indicated  (Reference  9).  This  strain  corresponds 
to  that  predicted  at  159  meters  for  the  Shoal  event. 

Additionally,  an  examination  of  uniaxial  compressive 
tests  on  Unaweep  granite  cores  (Reference  10)  revealed 
that  a  strain  of  0.33  percent  was  sufficient  to  produce 
fracturing.  Utilizing  this  strain,  a  free-field  fracture 
radius  of  146  meters  was  predicted  for  the  Shoal  event. 
These  predictions  appeared  reasonable  in  the  light  of  maxi 
mum  and  minimum  figures  previously  established  and  the 
fact  that  fracture  strains  for  various  granites  vary  only 
slightly.  A  radius  of  159  meters  was  therefore  predicted 
as  the  outer  limit  of  free-field  fracturing. 

2.2.3  Radius  of  Radiation  Imection  into  Cracks  by 
Particulates  and  Condensates.  Data  concerning  the  radial 
distribution  of  radioactive  materials  in  fractures  are 
very  limited.  The  outer  limit  of  radioactivity  indicated 
at  Rainier  was  about  5S  meters  (Reference  11).  The  outer 
limit  of  fracturing  shown  in  Table  14  (Reference  12)  was 
85.3  meters.  Assuming  that  the  ratio  of  the  radius  of 
radioactivity  to  the  radius  of  cracking  is  constant  for 
various  media,  then  the  outer  limit  of  radioactivity  for 
Shoal  would  have  been 

(58/85.3)  X  159  meters  =  108  meters 
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where  159  meters  is  the  cracking  radius  as  predicted  in 
Section  2.2.2.  A  check  was  made  by  considering  that  ra¬ 
dioactivity  entered  the  Hardhat  drift  at  a  distance  of 
at  least  61  meters.  By  a  computation  similar  to  that 
outlined  in  Section  2.2.2  a  minimum  radioactive  fracture 
radius  for  Shoal  was  found  to  be  about  91.4  meters. 

2.2.4  Height  of  Chimney.  A  chimney  height  of  110 
meters  was  predicted.  This  was  based  on  the  calculation 
described  in  Appendix  A  for  a  predicted  cavity  radius  of 
24.7  meters.  The  swell  factor,  n,  used  in  the  computa¬ 
tion,  was  based  on  the  observed  chimney  volume  for  Haufd- 
hat  in  relation  to  the  in  situ  volume  of  rock  which  filled 
the  available  space  following  cavity  collapse.  This 
factor  was  1.35,  as  computed  from  the  equations  given 

in  Appendix  A. 

2.2.5  Depth  of  Spalling.  Predictions  of  the  lowest 
limit  of  spalling  for  the  Shoal  event  were  obtained  by 
analog  computation.  The  spall  depths  and  thickness  of 
the  relatively  undisturbed  section  predicted  therefore 
differed  slightly  from  that  previously  quoted  (Reference 
13)  but  the  conclusions  as  to  containment  remained  valid. 

Since  the  number  and  depth  of  aaiy  spalls  calculated 
are  extremely  dependent  on  the  tensile  strength  that  is 
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assumed  for  the  medium,  various  cases  wert'  analyzed  by 
analog  computations  for  the  Shoal  event.  A  peak  pressure 
pulse  reflected  at  the  surface  equal  to  210  bars  (re¬ 
sulting  from  a  yield  of  14,26  kt )  was  assumed,  and  the 
tensile  strength  was  varied  from  6.9  to  82.7  bars  in 
6.9  bar  steps  to  determine  its  effect  on  spalling  depth. 

The  peak  pressure  due  to  a  12,5  kt  yield  varies  only 
slightly  from  the  above,  so  the  cal  cul  a  ti  ons  i<'inained  valid. 

The  graphical  plot  of  the  data  (Figure  2,3)  shows 
the  variation  in  depth  and  number  of  spalls  as  a  function 
of  tensile  strength.  It  may  be  seen  that  no  spalls  occur 
when  the  tensile  strength  exceeds  approximately  75.9  bars, 
and  the  deepest  spalling  occurs  when  the  assumed  strength 
is  around  13.8  bars.  The  analog  computations  were  checked 
and  found  to  agree  with  those  obtained  graphically  for 
27,6  bars.  One  could  then  have  said  that,  regardless  of 
the  tensile  strength,  the  worst  depth  of  spalling  should 
be  about  122  meters  below  the  surface.  This  lower  limit 
would  leave  a  relatively  undisturbed  section  of  134  meters 
between  it  and  the  top  of  the  predicted  chimney. 

2.2.6  Radius  of  Spalling.  To  determine  the  distance 
out  to  which  surface  damage  would  occur,  the  lateral  ex¬ 
tent  of  spalling  from  surface  zero  was  investigated.  The 
exact  determination  of  this  radius  was  not  obtainable  by 
analytic  methods,  but  estimates  were  obtained  by  comparing 
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the  predicted  surface  displacements  with  those  for  an 
event  where  the  spalling  radius  was  known .  For  Hardhat 
the  observed  surface  zero  displacement  due  to  spalling 
was  69.9  cm  and  spalling  extended  laterally  between  457 
and  914  meters  (Reference  14) .  Displacements  dropped 
off  rapidly  as  the  horizontal  distance  from  surface  zero 
increased.  At  457  meters  only  2.5  cm  was  observed.  For 
Shoal  the  surface  zero  displacement  was  predicted  to  be 
47.7  cm,  which  was  less  than  that  for  Hardhat.  Since  the 
decrease  in  surface  displacement  with  distance  was  as¬ 
sumed  to  be  similar  to  that  for  Hardhat,  the  lateral  ex¬ 
tent  of  spalling  was  predicted  to  be  of  about  the  same 
magnitude,  or  less  than  914  meters. 

2.2.7  Elastic  Radius.  The  elastic  radius  was  pre¬ 
dicted  to  occur  between  a  radius  of  253  to  286  meters. 
Prediction  of  the  elastic  radius  becomes  important  in  the 
determination  of  the  attenuation  rates  for  the  propagated 
seismic  signal.  From  an  examination  of  empirical  data, 
a  distinct  decrease  in  the  attenuation  is  noticed  on  ac¬ 
celeration  and  particle  velocity  records  at  a  particular 
distance  from  ground  zero.  This  point  may  be  accepted 
as  the  elastic  limit  and  may  be  determined  by  considering 
the  predicted  pressures  or  particle  velocities  in  relation 
to  the  properties  of  the  medium. 

Hydrodynamic  calculations  indicated  that  by  a  process 
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of  convergence  of  the  propagation  velocity  upon  a  final 
value  (i.e.,  the  limit  of  hydrodynamic  action),  an  clastic 
radius  would  be  determined  to  be  about  253  meters.  A 
similar  computation  for  Hardhat  yielded  a  radius  of  230 
meters,  which  was  in  general  agreement  with  that  observed. 
From  an  examination  of  the  regression  lines  for  particle 
velocity  at  Hardhat  by  the  Sandia  Corporation  (Reference 
15),  a  break  point  was  indicated  at  a  particle  velocity 
of  about  2.13  m/sec  and  a  range  of  222  meters.  A  possible 
additional  check  was  made  by  scaling  up  this  distance  to 
the  Shoal  yield  giving  a  prediction  of  286  meters  as  the 
elastic  limit. 

Some  doubt  may  be  entertained  as  to  whether  the 
clastic  limit  is  a  distinct  point  or  rather  a  range  of 
pressures  or  particle  velocities  over  which  a  change  in 
attenuation  occurs.  The  assignment  of  a  discrete  radius 
was  therefore  extremely  tenuous  and  the  elastic  limit 
was  considered  as  occurring  approximately  within  a  zone 
given  by  the  radii  shown  above  (253  to  286  meters). 

2.2.8  Tunnel  Closure  Radius.  Tunnel  closure  for 
the  Hardhat  event  was  continuous  from  the  end  of  the 
drift  out  to  station  3  90  (Reference  8),  a  radial  dis¬ 

tance  of  128  meters  from  the  W.P.  Using  this  value  as 
the  parameter  for  cube-root  scaling,  a  prediction  of  164 
meters  was  obtained  as  the  outer  limit  of  complete  tunnel 
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closure  for  Shoal.  As  in  the  case  of  continuous  tunnel 


damage,  closure  is  dependent  upon  the  pre-shot  local 
geology  to  a  large  extent.  Since  the  ground  out  to 
station  3+85  was  highly  fractured,  the  possibility  of 
closure  to  187  meters  was  not  considered  unlikely.  With 
this  in  mind,  and  the  fact  that  some  closure  at  Hardhat 
occurred  at  137  meters  in  a  fault  zone  (Reference  8), 
closure  for  Shoal  was  anticipated  to  a  radial  distance 
of  187  meters. 

2,2.9  Tunnel  Damage  Radius.  The  extent  of  damage 
to  tunnels  and  shafts  arising  from  a  subsurface  detona¬ 
tion  is  a  function  of  the  pre-shot  geologic  condition, 
physical  properties  of  the  mediiun,  the  angle  of  inci¬ 
dence  and  any  support  features  employed.  Radii  quoted 
in  Section  2.2.8  for  this  terminal  effect  considered 
these  factors  and  indicated  a  range  over  which  this 
damage  would  occur.  However,  since  it  is  more  useful 
to  indicate  a  specific  distance  and  the  factors  affecting 
possible  variations,  a  more  comprehensive  procetiure  was 
employed.  The  availability  of  empirical  data  from  Hard¬ 
hat,  an  event  in  granite  with  an  essentially  similar 
configuration,  provided  the  basis  for  orediction  of 
effects  bv  cube-root  scaling  witnout  tne  omission  of 
extenuating  j.ac\o 

Continuous  damage  in  the  Hardnat  access  tunnel 
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extended  to  station  2  +  47,  a  radial  distance  of  about 
172  meters  from  the  shot  point  (Reference  8) .  The  heavy 
damage  between  stations  2  ■♦•'‘90  and  3  40  was  attributed 

to  incompetent  pre-shot  conditions  produced  by  a  fault 
zone  in  this  region.  The  radius  of  continuous  tunnel 
damage  for  Hardhat  was  172  meters.  Using  this  value, 
by  cube-root  scaling,  a  damage  distance  of  221  meters 
was  predicted  for  Shoal.  It  was  considered  possible  that 
the  presence  of  a  highly  fractured  and  faulted  zone  at 
a  range  of  232  to  258  meters  might  increase  this  predicted 
distance  but  the  amount  was  not  subject  to  accurate  com¬ 
putation  . 

2  2.10  Predictions  on  Stemming.  The  evaluation  of 
the  stemming  plan  for  fhoal  contained  in  this  report  was 
based  on  the  information  contained  in  Reference  16.  Con¬ 
sideration  was  also  given  to  further  details  contained 
in  other  construction  and  support  criteria.  The  yield 
used  was  the  measured  value  of  12.5  kt  as  stated  in 
Reference  17.  The  depth  of  W.P.  was  366  meters,  and 
the  location  was  6.1  meters  south  of  core  hole  ECH~D, 

At  the  NVOO  meeting  on  acceptability  September  7,  1962, 
the  criterion  established  by  NVOO  was  that  venting  should 
not  exceed  that  of  the  Gnome  and  Hardhat  events. 

Since  the  Hardhat  event  was  the  only  IJ .  S .  nuclear 
explosion  in  granite,  a  comparison  was  made  the 
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effects  of  that  event  and  those  which  were  anticipated 
for  the  Shoal  event.  In  both  cases  the  chamber  sizes 
were  very  small.  Accordingly,  the  effect  of  chamber 
volume  was  not  an  important  factor  in  comparing  the  two 
events . 

The  design  for  emplacement  of  the  device  vertically 
over  the  access  drift  instead  of  to  one  side  was  new. 

Also,  there  were  a  number  of  features  of  particular  in¬ 
terest.  The  device  was  placed  in  a  38.1  cm  diameter  hori¬ 
zontal  hole  6.1  meters  from  the  W.P.  and  stemmed  with 
sand.  This  arrangement  was  considered  comparable  to  the 
explosion  of  a  device  in  a  drilled  hole.  The  stemmed 
hole  and  sand  near  the  W.P.  was  considered  to  provide 
effective  initial  stemming. 

The  back  or  top  of  the  access  tunnel  directly  be¬ 
low  the  proposed  device  emplacement  locations,  Station 
10  +  00,  was  at  a  distance  of  9.15  meters  from  the  center 
of  the  explosion.  Calculations  indicated  that  the  tunnel 
would  be  closed  before  any  blast  effects  could  reach  this 
location  by  the  longer  path  through  Plug  #1 ,  down  the 
raise  and  back  to  Station  10  OO.  The  fact  that  sand 
stemming  originally  proposed  in  the  raise  was  omitted 
was  not  expected  to  adversely  affect  the  performance  of 
the  self-sealing  arrangement. 

The  minimum  yield  that  would  produce  closure  was 
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oosnpuit®xi  on  the  basis  of  equal  Impulse  to  be  0.02  kt.  The 
usethods  of  computations  for  such  small  yields  were  open 
to  cuestlon.  However,  even  if  the  yield  were  very  small, 
it  was  considered  that  the  series  of  plugs  would  be  ade¬ 
quate  to  contain  the  explosion. 

2,2„11  Effectiveness  of  Plugs  in  Access  Drift.  The 
eff ectiveness  of  plugs  is  related  to  the  damage  radii  an¬ 
ticipated.  The  four  sand-fill  plugs  were  at  distances 
along  the  di'ift  of  approximately  6.1,  27.4,  45.7  and  76.2 
meters  from  the  point  opposite  the  W.P.  The  latter  two 
were  much  closer  scaled  radial  distances  than  the  sand 
plugs  at  Hardhat.  These  were  located  at  approximate 
distances  of  84  and  110  meters.  Experience  at  Hardhat 
relative  to  the  latter  plug  Indicated  that  inward  move¬ 
ment  of  rock  compressed  the  plug  thereby  increasinr  it* 
effectiveness.  However,  following  cavity  collapse  j 

predicted  for  Shoal),  radioactive  gases  permeated  the 
broken  granite  and  entered  the  access  tunnel  beyond  the 
last  plug.  These  gases  were  largely  contained  by  the 
gas  seal  door  which  was  located  approximately  234  meters 
from  the  W.P.  At  Shoal  similar  effects  were  considered 
probable  but  Plug  #7  in  the  shaft  was  felt  to  be  suffi¬ 
cient  to  filter  out  any  radioactive  particulate  matter 
that  might  reach  this  area  and  thus  minimize  escape  of 
radioactive  gases  to  the  surface. 
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2.2.12  Conclusions.  As  a  result  of  evaluation  of 


Reference  16  and  subsequent  data,  conclusions  were  as 
follows ; 

1.  There  would  be  damage  to  the  access  drift  beyond 
the  last  plug  in  the  access  drift  with  a  high 
probability  of  contamination  of  that  portion  of 
the  drift . 

2.  The  30.5  meter  sand  plug  at  the  bottom  of  the 
shaft  would  prevent  radioactive  fiarticulate 
fnatter  from  reaching  that  portion  of  the  shaft 
above  the  plug  and  reduce  the  venting  of  radio¬ 
active  gases  to  an  acceptable  level. 

3.  Damage  to  the  shaft  would  necessitate  remedial 
measures  comparable  to  those  for  Hardhat  if 
re-entry  were  made. 

4.  Geologic  conditions  were  less  favorable  for  con¬ 
tainment  than  at  Hardhat  in  that  faulting  and 
fracturing  were  extensive.  This  would  be  a 
factor  in  damage  generally  increasing  damage 
where  faults  exist. 

5.  The  stemming  plan  was  generally  adequate  to  re¬ 
duce  venting  to  a  level  less  than  observed  at 
the  Gnome  and  Hardhat  events. 

2.2.13  Evaluation  of  Possibility  of  St  earn -Generated 
Secondary  Explosion.  Since  hydrologic  conditions  at  Shoal 
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were  not  accurately  knowm  (see  page  14) ,  thp  presence  of 
water-bearing  zones  in  the  vicinity  of  the  W.P.  was  as¬ 
sumed  and  possible  effects  were  investigated.  A  review 
of  the  probable  temperature-pressure  history  of  the  cavity 
was  first  made.  This  showed  that  initially  the  cavity 
would  be  filled  with  vapx>rized  granite  and  water  in  the 
form  of  steam,  the  water  being  obtained  from  fixed  water 
in  the  granite,  all  at  extremely  high  temperature  and 
pressure.  Almost  Immediately  the  temperature  would  Irop 
to  a  level  where  vaporized  granite  would  condense.  At 
the  end  of  this  phase  the  cavity  would  be  filled  with 
steam  and  there  would  be  molten  granite  on  the  side  walls 
and  probably  collected  in  a  pool  at  the  bottom.  The 
temperature  would  then  continue  to  be  reduced  by  heat 
transfer  to  surrounding  rock  until  a  plateau  was  reached, 
when  the  granite  started  to  solidify.  The  reason  for 
the  plateau  is  that  as  the  granite  solidifies,  it  gives 
up  its  heat  of  fusion,  slowing  down  the  general  cooling. 
This  temperature  would  be  about  704  degrees  C,  the  ap¬ 
proximate  melting  point  of  granite  (Reference  18). 

Consideration  was  then  given  to  the  situation  ex¬ 
isting  when  the  temperature  remains  at  704  degrees  C 
during  the  period  in  which  solidification  of  molten 
granite  was  procecKllng  at  a  relatively  slow  rate.  As 
pointed  out  above,  all  the  vaporized  rock  (and  metal 
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from  the  device)  would  have  been  condensed.  The  device 
was  placed  in  a  38.1  cm  hole  which  was  tamped;  therefore, 
there  should  be  very  little  air  in  the  cavity.  Only  steam 
produced  from  water  contained  in  vaporized  or  molten  rock 
would  remain.  It  was  believed  that  the  fixed  water  in  the 
rock  which  was  displaced  radially  would  not  be  released. 
Assuming  that  cavity  collapse  had  not  occurred,  a 
computation  was  made  as  shown  in  Appendix  B.  This  was 
for  a  25.9  meter  radius  cavity  but  was  also  valid  for  a 
24.7  meter  radius  predicted  for  a  yield  of  12.5  kt.  Of 
particular  interest  was  the  conclusion  that  after  vapor¬ 
ized  rock  has  condensed,  the  cavity  is  filled  with  steam 
at  704  degrees  C  and  the  pressure  would  only  be  .345  bar 
absolute.  In  Reference  19,  a  post-shot  evaluation  of 
Rainier,  it  is  hypothesized  that  a  pressure  drop  occurred 
because  steam  escaped  from  the  cavity  through  fissures. 
(Since  there  was  no  report  of  venting,  this  assumption 
of  escape  of  steam  may  not  be  justified.)  The  occurrence 
of  the  pressure  drop  is  based  on  examination  of  vesicu- 
lated  material.  It  is  stated  on  page  4,  Reference  19, 
that,  ’’Thus,  the  evidence  strongly  points  to  the  fact 
that  the  vapor  pressure  was  40  atmospheres  inside  the 
chamber  immediately  after  the  shot  while  glass  droplets 
were  still  falling  from  the  roof  of  the  chamber;  a  minute 
or  two  later  the  pressure  had  dropped  to  a  few  atmospheres." 
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The  drop  of  pressure  to  a  low  valxie  immediately  after  the 
explosion  was  considered  to  be  a  factor  in  cavity  collapse. 

The  spalling  of  pieces  of  rock  from  the  surface  of 
the  cavity  by  steam  forming  behind  the  face  would  be 
accelerated  by  the  low  pressure.  However,  the  analysis 
indicated  that  an  influx  of  a  large  volume  of  water  with¬ 
in  a  particular  length  of  time  would  be  required  to  cause 
the  pressure  to  rise  to  6.9  bars.  At  the  assumed  rate 
of  flow  of  265.0  liter/min,  7  hours  would  be  required. 

If  a  larger  flow  was  assumed  the  time  would  be  reduced 
proportionately.  Although  precise  procfedu.res  for  analy¬ 
sis  of  the  rates  of  heat  transfer  were  not  developed,  an 
approximate  analysis  indicated  that  the  inflow  of  water 
would  have  a  net  cooling  effect  on  the  heated  rock  in 
and  surrounding  the  cavity. 

The  above  discussion  assumed  that  the  cavity  would 
not  collapse  before  the  pressure  was  reduced  to  the  point 
where  water  might  enter.  Since  at  Hardhat ,  cavity  col¬ 
lapse  did  not  occur  until  11  hours  after  the  event  (Ref¬ 
erence  7),  this  was  a  reasonable  assumption.  It  was 
assumed  that  cavity  collapse  was  unlikely  until  the  pres¬ 
sure  in  the  cavity  was  reduced  to  a  value  less  than  the 
lithostatic  pressure,  82.7  bars.  If  collapse  occurred 
when  the  cavity  pressure  was  somewhat  greater  than  6.9 
bars,  cooling  should  have  proceeded  at  a  rapid  rate  from 
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the  effect  of  the  relatively  cool  rock  masses  falling  in¬ 
to  the  cavity. 

The  only  other  concern  was  that  water  entering  tun¬ 
nels  in  the  vicinity  of  the  shot  point  may  have  been 
heated  by  the  effect  of  the  explosion.  Consideration  of 
the  limited  heat  available  in  such  a  case  indicated  that 
relatively  small  quantities  of  steam  would  be  generated. 
While  such  an  explosion  might  have  had  a  local  effect, 
it  was  felt  unlikely  that  it  could  affect  containment. 
Passage  of  the  steam  through  access  tunnel  plugs  was  con¬ 
sidered  improbable.  If  it  did  vent  into  the  access  tun¬ 
nel,  condensation  would  reduce  the  pressure  materially 
before  the  steam  reached  the  30.5  meter  plug  in  the  shaft. 
This  massive  plug  was  felt  sufficient  to  effectively  pre¬ 
vent  any  steam  escaping  to  the  atmosphere. 

On  the  basis  of  the  above  discussion,  the  following 
conclusions  were  reached: 

1.  Evidence  suggested  that  ground  water  may  have 
existed  in  the  vicinity  of  the  W.P. 

2.  Following  the  condensation  of  vaporized  rock 
and  metal  from  the  device,  the  steam  pressure 
in  the  cavity  was  expected  to  be  very  low,  less 
than  one  bar.  When  the  shell  of  melted  rock 
was  broken  by  Implosion  of  portions  of  the  shell 
by  steam  in  pockets  behind  the  face,  an  inflow 
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of  water  would  occur  slowly,  raising  the  steam 
pressure  and  cooling  the  rock.  The  pressure  was 
considered  to  be  limited  to  6.9  bars  as  any  in¬ 
crease  above  that  pressure  would  prevent  further 
inflow  of  water. 

3.  Shourd  the  cavity  collapse,  temperatures  and 
pressures  would  be  lower  than  at  the  time  rock 
is  solidifying.  Rapid  cooling  from  introduction 
of  broken  rock  would  reduce  temperatures  and 
pressures  rapidly, 

4.  Entry  of  water  into  heated  tunnels  or  other 
openings  near  the  W.P.  might  cause  localized  ex¬ 
plosions.  Condensation  of  steam  in  passageways 
and  the  30.5  meter  plug  in  the  shaft  would  have 
prevented  steam  reaching  the  surface. 

5.  Finally,  no  significant  hazard  existed  from 
secondary  explosion  by  entry  of  ground  water 
into  the  cavity. 

2,2.14  Estimate  of  Prompt  Venting  of  Gaseous  Isotopes. 
It  was  estimated  that  prompt  venting  of  radioactive  iso¬ 
topes  to  the  surface  through  fissures  was  unlikely  for 
reasons  stated  below.  The  possibility  of  venting  through 
stemming  is  discussed  in  Section  2.2.12. 

It  was  predicted  that  surface  spalling  would  extend 
to  a  depth  at  Surface  Zero  not  greater  than  122  meters. 
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Prior  to  the  collapse  of  the  cavity  (which,  on  the  basis 
of  experience  in  Hardhat,  would  be  delayed)  there  would 
remain  above  the  cavity  a  thickness  of  366  -  (25.9  +  122)  = 
218  meters  of  relatively  undisturbed  rock.  All  holes  from 
the  surface  were  reported  to  have  been  grouted  to  prevent 
escape  of  radioactive  gases  to  the  surface.  Drilling  in 
the  vicinity  of  the  W.  P.  established  the  fact  there  were 
no  faults  within  a  radial  distance  of  30.5  meters.  Under 
these  conditions  the  displacement  of  material  radially 
from  the  W.  P.  would  compress  the  gouge  in  the  fault  and 
inhibit  the  escape  of  radioactive  gases  through  Fault 
Zone  E  at  a  radial  distance  of  about  34  meters. 

The  possibility  of  immediate  escape  of  radioactive 
gases  through  a  thickness  of  218  meters  appeared  to  be 
extremely  remote.  In  the  first  place  the  very  small  di¬ 
ameter  of  the  chamber  assured  the  maximum  degree  of  melting 
of  the  granite.  The  molten  layer  at  the  periphery  of  the 
cavity  was  felt  to  be  Impervious  to  the  escape  of  gases 
as  the  tendency  would  be  to  drive  the  molten  material  in¬ 
to  any  fissures  which  might  exist  and  close  them.  Such 
action  is  indicated  in  tuff  (Reference  19).  However,  as 
cooling  occurs,  post -shot  observations  of  Rainier  indicate 
that  the  wall  material  is  Imploded  into  the  cavity  (Ref¬ 
erence  19).  This  is  believed  to  be  caused  by  steam  pres¬ 
sure  behind  the  surface.  Since  granite  has  a  much  smaller 
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quantity  of  moisture  and  greater  competence,  it  was  con¬ 
sidered  unlikely  that  such  implosion  would  occur  to  the 
same  extent  or  as  rapidly.  This  is  substantiated  by  the 
fact  that  no  radiation  was  detected  at  Rainier 's  Surface 
Zero  until  H  1  hour  at  which  time  it  was  only  lO  mr/hr 
(Reference  19). 

As  a  further  consideration,  estimates  of  the  small 
amount  of  water  available  for  generation  of  steam  indi¬ 
cated  that  pressures  Ir  the  cavity  should  be  very  low 
following  condensation  of  gaseous  rock,  which  should 
occur  rapidly.  This  was  Important  as  there  would  be 
insufficient  pressure  in  the  cavity  to  force  radioactive 
gases  up  through  the  overlying  rock  (see  Section  2.2,13), 

2,2,15  Prediction  of  Containment.  It  was  predicted 
that  the  Shoal  event  would  be  contained  to  the  same  degree 
as  the  Hardhat  event  for  reasons  discussed  below. 

Containment  of  the  explosion  at  the  time  of  the 
event  depends  upon  adequacy  of  stemming,  depth  of  burst, 
and  absence  of  passages  which  could  provide  a  path  for 
radioactive  venting.  Following  the  explosion,  some  of 
the  wall  material  was  expected  to  implode,  br*'aking  the 
seal  formed  by  the  melting  of  a  shell  of  roeV  at  the 
periphery  of  the  cavity.  This  would  permit  i  lio  escape 
of  some  radioactive  gases  into  fissures  in  tlu  rock. 

After  the  cavity  collapsed,  a  different  ■situation 
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would  exist.  A  chimney  having  a  height  '*<  I  If)  meters  mea¬ 
sured  from  the  elevation  of  the  W.P.  was  tjredicteci.  A 
depth  of  spalling  not  to  exceed  122  meters  was  predicted. 
Thus,  there  would  remain  a  thickness  of  3f>6  -  (110  +  122) 

=  134  meters  of  relatively  undisturbed  rock  as  a  barrier 
against  escape  of  radioactive  gases  to  the  surface.  It 
was  assumed  that  the  temperatures  and  pressures  would  have 
been  markedly  reduced  by  the  collapse  of  the  cavity  and 
the  filling  of  the  cavity  with  rubble  at  the  ambient  tem¬ 
perature  of  the  rock  (except  for  that  portion  adjacent  to 
the  cavity).  Therefore,  any  possible  escape  of  radioactive 
gases  to  the  surface  would  be  more  in  the  nature  of  a  low- 
pressure  migration  or  diffusion.  Radioactive  isotopes 
that  condense  at  ambient  temperatures  should  not  reach 
the  surface.  There  remain  xenon  and  krypton  which  might 
possibly  escape  to  the  surface  at  a  slow  rate. 

Stemming,  as  discussed  in  Section  2.2.10,  was  pre¬ 
dicted  to  be  adequate  for  the  immediate  period  of  the 
explosion  and  that  following  cavity  collapse. 

2.3  SEISMIC  EFFECTS 

Predictions  of  particle  motions,  for  the  Shoal  model 
described  in  Chapter  1,  were  made  using  the  equations 
given  below. 

2.3.1  Acceleration.  The  particle  acceleration 
equation  was  developed  through  use  of  a  1  east-squares 
equation  for  three  components  of  surface  peak  particle 
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acceleration  for  tuff  which  was  adjusted  to  fit  the  value 
of  acceleration  observed  at  the  Las  Vegas  Seismic  Safety 
Net  Station  from  the  Clearwater  event. 

The  least-squares  equation  was: 

a  -  3.2  X  10^  R"^*®  (Reference  20)  (2.3.1) 

where 

a  is  in  units  of  gravity 
W  is  in  kilotons 
R  is  in  feet 
or  for  R  in  meters, 

a  =  2.96  X  10^  w'^  r“^‘®  (2.3.2) 

Computation  of  peak  surface  particle  acceleration 
using  equation  2.3.2,  the  Clearwater -Las  Vegas  distance 
in  meters  and  the  reported  yield  for  Clearwater  yielded 
a  value  lower  than  the  recorded  0.0005  g  (Reference  21). 
To  rectify  the  discrepancy,  the  coefficient  of  equation 
2.3.2  was  adjusted  to  fit  the  observed  value  by  sub¬ 
stituting  the  observed  Clearwater -Las  Vegas  results  in 
the  equation: 

a  =  K  R"^*^  (2.3.3) 

The  value  of  the  coefficient,  K,  obtained  from  equa¬ 
tion  2.3.3  is  5.03  x  lO^  which,  when  substituted  above, 
yields  the  equation: 

a  =  5.03  X  10^  R“^'®  (2.3.4) 

This  was  used  to  predict  surface  particle  acceleration 
for  the  Shoal  event  (see  Figure  2.4  and  Table  2.2). 
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2,3,2  Displacement. 


The  displacement  equation  is 
that  reported  in  Reference  22  as: 

Log  A  =  0.75  Log  Y  -  1.73  Log  R  (2.3.5) 

which,  rewritten  in  exponential  form,  is: 

d  =  (2.3.6) 

where 

d  =  surface  displacement  in  centimeters 

W  =  yield  in  kilotons 

R  =  distance  in  kilometers 
or  for  R  in  meters, 

5  75-1.73  .  . 

d  =  1.51  X  10^  R  (2.3.7) 

The  peak  surface  particle  displacement  predicted  at 
Las  Vegas  for  the  Clearwater  event  was  0.004  cm.  Since 
a  value  of  0.003  cm  (Reference  21)  was  observed  for  a 
lower  yield,  equation  2.3.7  was  used  to  predict  the  dis¬ 
placements  expected  from  the  Shoal  event  (see  Figure  2.5 
and  Table  2,2) . 

2.3,3  Velocity.  Data  for  peak  surface  particle 
velocity,  of  a  reliability  equal  to  that  of  acceleration 
and  displacement  data  in  the  range  of  interest,  were  not 
available.  Accordingly,  simple  harmonic  motion  was  as¬ 
sumed  and  an  equation  for  particle  velocity  was  obtained 
through  use  of  the  relationship} 


where 
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K  =  gravitational  conotant  ~  980  c^nt iinet  ers/second/ 
second 

u  -  peak  surface  particle  /elocity  in  contimeters/ 
second 


a  =  peak  surface  particle  acceleration  in  g’s 
d  =  peak  surface  p>article  diisplacement  in  centimeters. 
The  values  of  acceleration  and  displacement  substi¬ 
tuted  in  equation  2.3.8  are  obtained  from  prediction  equa¬ 
tions  2.3.4  and  2.3.7  given  above.  The  resulting  peak 
surface  particle  velocity  equation  is: 

u  -  8.64  X  10®  R  (2.3.9) 

This  equation  was  used  to  predict  particle  velocity 
for  Shoal  (see  Figure  2.6  and  Table  2.2). 


2.3.4 


mat ions  in  Terms 


Particle  Motion.  Equa¬ 


tions  2.3.4,  2.3.7  and  2.3.9  may  be  rewritten  in  order  to 


obtain  distance  to  any  particular  motion  of  interest.  These 


are: 


j^5.03  X  10^  j  ' 
jl.51  X  10^  * 
je.64  X  10^  ' 


(2.3.10) 

(2.3.11) 

(2.3.12) 


where 

R  s  distance  in  meters 
W  s=  yield  in  kilotons 
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a  =  peak  surface  particle  acceleration  in  units  of 
gravity  (g) 

d  =  peak  surface  particle  displacement  in  centimeters 
u  =  peak  surface  particle  velocity  in  centimeters/ 
second 

Equations  2,3.10  and  2.3.12  have  been  used  in  com¬ 
puting  distance  to  0.1  g  and  11  ctn/sec  for  radius  of  damage 
as  shown  in  Figure  2.7  and  distance  to  0.001  g  for  limit 
of  perceptibility  in  the  same  figure. 
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CHAPTEH  1 1 1 


ANALYSIS  AND  INTERPRETATION 

3.1  PHENOMENOLXWY  AND  CONTAINMENT 

As  stcited  in  Chapter  2,  prediction  of  effects  and  an 
evaluation  of  containment  were  made  for  the  purposes  of 
operational  safety.  Limited  post-shot  exploration  permits 
a  comparison  of  only  some  of  the  predicted  quantities  with 
observed  values.  The  primary  result  of  the  Shoal  event 
was  that  it  was  satisfactorily  contained  as  predicted. 

The  prediction  of  containment  was  based  on  an  estimation 
and  analysis  of  the  individual  effects  anticipated  and  the 
interrelation  of  one  to  another  with  respect  to  containment. 

3.1.1  Cavity  Radius.  The  post-shot  cavity  radius  as 
determined  by  drilling  was  25,6  +  .6  meters  (Reference  17). 
This  compares  very  well  with  the  predicted  radius  of  24.7 
meters,  confirming  the  procedure  of  scaling  Hardhat  results 
to  the  Shoal  yield. 

3.1.2  Chimney  Height.  The  post-shot  chimney  height 
was  found  to  be  100  meters  (Reference  17).  The  predicted 
value  was  110  meters.  This  close  agreement  substantiates 
the  granite  swell  factor  and  other  assumptions  used  to  de¬ 
termine  the  extent  of  chimney  formation. 
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3.1.3  Depth  of  Spalling.  The  lowest  limit  of  spall- 
tig  was  predicted  to  be  122  meters  below  the  surface.  Ex¬ 
act  post-shot  data  are  not  available  but  indications  arc 
that  spalling  extended  below  120  meters.  It  has  been 
stated  (Reference  23)  that  an  examination  of  records  from 
gages  vertically  over  the  shot  point  indicated  that  free 
fall  occurcd  at  a  depth  of  85.5  meters.  It  was  also  be¬ 
lieved  that  free  fall  occurred  at  a  depth  of  120  meters 
but  owing  to  gage  difficulties  this  is  not  positive.  The 
cable  leading  to  the  accelerometer  broke  at  0.3  seconds 
providing  only  a  partial  record.  The  velocity  gage  gave 
a  continuous  record  and  the  slope  indicated  free  fall  but 
the  period  was  so  short  that  a  positive  interpretation  can- 
no  t  be  made . 

3.1.4  Radius  of  Spalling.  The  lateral  extent  of 
spalling  from  surface  zero  was  predicted  to  be  less  than 
914  meters.  This  has  been  shown  to  be  close  to  the  ob¬ 
served  radius  by  the  results  given  in  Reference  24.  From 
aerial  examination  it  was  found  that  intense  fracturing 
was  confined  to  within  107  meters  of  surface  zero.  Frac¬ 
tures  up  to  1.3  cm  wide  occurred  out  to  458  meters  south, 
794  meters  north,  366  laeters  east  and  610  meters  west; 
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some  fractures  in  this  zone  were  as  wide  as  15  cm.  It  is 


felt  that  these  fractures  are  indicative  of  the  radii  to 
which  spalling  occurred  and  are  in  agreement  with  the  pre¬ 
diction  thus  providing  substantiation  for  the  procedures 
used . 

The  above-mentioned  effects  were  all  used  to  arrive 
at  the  conclusion  that  the  Shoal  event  would  be  contained. 
Additionally,  an  evaluation  of  the  stemming  plan  showed  it 
to  be  adequate.  Both  of  these  conclusions  were  completely 
substantiated  by  the  results  of  the  test. 

Several  predicted  effects  such  as  cracking  radius, 
radius  of  radioactive  injection  into  cracks  and  tunnel 
damage  predictions  c^^n  not  be  substantiated  accurately 
due  to  the  lack  of  post-shot  data  at  this  time. 

3.2  SEISMIC  EFFECTS 

As  detailed  in  Chapter  1,  surface  motions  resulting 
from  the  Shoal  event  were  recorded  at  various  locations 
with  acceleration,  displacement,  and/or  velocity-sensitive 
instruments.  The  resulting  records  have  been  corrected 
for  the  particular  Instrument  transfer  function  to  obtain 
true  ground  motion  neglecting  instrument  to  ground  cou¬ 
pling.  Additionally,  the  corrected  acceleration  and 
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displacement  records  have  been  integrated  and  differen¬ 
tiated  respectively  to  obtain  surface  particle  velocity 
where  direct  measurements  of  this  parameter  were  not 
made.  Peak  values  of  surface  particle  motion  have  been 
determined  from  the  corrected  records.  In  addition,  the 
available  components  at  each  station  have  been  analyzed 
simultaneously  in  order  to  determine  the  absolute  values 
of  peak  motion  received  at  the  station.  These  values  are 
denoted  as  peak  resultant  vectors. 

In  order  to  evaluate  the  specific  objectives  stated 
in  Chapter  I ,  the  peak  motions  have  been  segregated  into 
classifications  on  the  basis  of  parameter  of  motion  meas¬ 
ured,  component  measured,  station  lithology  and  relation 
of  station  location  to  areal  geology.  All  peak  motions, 
both  directly  measured  and  derived,  are  given  in  Tables 
3.1  through  3.3.  In  view  of  the  fact  that  funds  for  ac¬ 
complishment  of  detailed  analyses  of  data  were  delayed, 
it  is  not  practicable  to  Include  such  results  in  this  re¬ 
port.  A  research  report  covering  the  analyses  of  Shoal 
data  and  that  from  other  events  will  be  submitted  when 
detailed  analyses  have  been  completed.  The  preliminary 
analysis  and  interpretation  of  results  presented  here 
are  given  on  an  objective-by-objective  basis. 
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3.2.1  The  Variation  in  Amplitude  of  Ground  Motions 


on  Granite  and  Alluvium  with  Distance.  Figures  3,1 
through  3.12  show  peak  surface  resultant  vectors  and  in¬ 
dividual  components  of  ground  motion  recorded  on  granite 
and  alluvium  from  the  Shoal  event  by  the  USC&GS.  Besides 
the  recorded  values,  the  graphs  shocv  the  least-square  re¬ 
gression  equation  derived  from  the  data  and  the  standard 
error  of  estimate  (which  will  be  discussed  in  detail  in 
Section  3.2.6,  below). 

Th^  data  recorded  on  aliuviua«  show  consistently 
larger  particle  motions  than  that  recorded  on  granite. 

The  rates  of  attenuation  for  the  two  media  are  nearly  the 
same.  Sufficient  measureicents  were  made  to  allovt:  an  ap- 
proximation  of  the  magnitude  of  amplification  caused  by 
the  presence  of  alluvium.  Although  soae  differences  in 
rates  of  attenuation  between  granite  and  alluvium  meas¬ 
urements  are  evident  (see  Figures  3.1  throagh  3.4)  and 
the  overlap  of  data  with  distance  is  small,  it  can  be 
seen  that  peak  accelerations  on  alJuviura  are  larger  than 
peak  accelerations  on  granite  by  an  average  factor  uf  2,.l. 
Similarly,  for  displacement  (Figures  3. .3  thr<iugh  3.8)  the 
factor  averaged  3.7  and  for  velocity  (Figures  3.9-3,12) 
3,4. 
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3.2.2  The  Effect  of  Faulting  on  Ground  Motion.  Peak 
resultant  vectors  of  particle  motion  toeasured  on  instru¬ 
ment  lines  parallel  to  the  major  northeast-southwest  fault 
trend  have  been  compared  with  those  measured  on  lines  per¬ 
pendicular  to  the  faulting.  The  orientation  of  applicable 
stations  with  respect  to  the  fault  trend  are  listed  in 
Table  3.4.  Figures  3.13  through  3.15  show  data  points  and 
the  associated  least-squares  regression  lines  for  lines 
parallel  and  perpendicular  to  the  fault  trend.  These  fig¬ 
ures  show  a  trend  to  higher  rates  of  attenuation  perpen¬ 
dicular  to  the  faulting  than  are  seen  parallel  to  the 
faulting.  However,  the  evidence  is  not  sufficiently 
strong  to  state  that  there  is  a  definite  difference  in 
attenuation  due  to  energy  transmission  across  fault  sys¬ 
tems  . 

3.2.3  The  Effect  of  Thickness  of  Alluv ium  on  the 
Amplitude  and  Attenuation  of  Amplitude  with  Distance. 

This  effect  may  be  best  examined  for  Shoal  by  comparing 
data  recorded  on  the  thick  alluvium  of  Fairview  Valley 
to  the  east  of  the  Sand  Springs  Range  with  that  recorded 
on  Fourmile  and  Eightmile  Flats  to  the  west  cf  the  Sand 
Springs  Range  where  the  alluvium  wedges  out  against  the 
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underlying  granite.  However,  insufficient  data  were  ob¬ 
tained  to  make  a  meaningful  statistical  evaluation  of 
this  effect. 

3,2.4  The  Partition  of  Compressional  Wave  Energy 
.  '<to  Surface  Waves  at  a  Lithologic  Discontinuity.  As 
een  in  Figures  3.1  through  3.12,  the  rates  of  attenua- 
tioii  for  stations  on  alluvium  are  nearly  the  same  as  those 
for  stations  on  granite.  This  seems  to  indicate  that 
even  though  thick  alluvium  is  present,  the  peak  motions 
measured  on  alluvium  were  the  result  of  the  same  type  of 
wave  that  produced  peak  motions  at  the  granite  stations. 
The  rates  of  attenuation  observed  on  alluvium  from  the 
Shoal  event  are  considerably  higher  (with  the  exception 
of  the  transverse  component)  than  the  attenuation  rate 
observed  for  surface-type  waves  generated  by  twenty-eight 
events  at  Yucca  Flat  (Reference  25).  Therefore,  it  is 
probable  that  all  peak  motions  recorded  at  Shoal  were  the 
result  of  body  waves  which  in  the  presence  of  thick  al¬ 
luvium  have  been  aoiplified  by  the  factors  discussed  in 
Section  3.2.1,  above.  Because  of  this,  no  determination 
of  energy  partition  between  body  and  surface  waves  at  a 
lithologic  discontinuity  can  be  made.  It  is  important 
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to  note  that  the  above  discussion  is  based  on  p>eak  arnpli 
tudes  measured  by  instruments  which  respond  best  to  peri 
ods  of  four  seconds  or  less.  The  intent  is  not  to  say 
that  long  period  surface  waves,  possibly  of  large  ampli¬ 
tude,  did  not  occur;  but  rather  that  there  was  no  evi¬ 
dence  from  these  particular  instruments  that  this  type 
wave  produced  any  peak  motions. 

3.2.5  The  Directional  Response  Characteristics  of 
the  Shot  Medium.  The  azimuthal  variations  of  amplitude 
in  the  Sand  Springs  Range  appear  to  be  controlled  to  a 
certain  extent  by  the  major  fault  trend  in  the  range  (see 
Section  3.2.2,  above).  Analysis  has  not  indicated  any 
azimuthal  variations  due  to  inherent  properties  of  the 
shot  media,  i.e.,  granite,  separately  from  the  effects 
attributed  to  faulting. 

3.2.6  The  Degree  of  Scatter  in  Recorded  Ground 
Mot  ions .  The  degree  of  scatter  is  difficult  to  judge 
without  the  use  of  statistical  analysis.  A  convenient 
numerical  quantity  is  the  standard  error  of  estimate 
(rfy)  which  can  be  determined  during  the  regression  an¬ 
alysis.  This  is  a  multiplicative  factor  --  additive  in 
logarithmic  form  --  and  has  the  following  significance. 
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Assuming  ntirmaJ  d  is  t  r  ibu  t  i  tm  of  data  points  usoci  in  the 
regrcsrion  .nalysis,  a  band  width  of  1  ^  ^  on  I’ither  side 

of  the  regression  line  will  include  about  b8  piercent  of 
the  data  points.  A  band  of  +  2  ay  or  3  rr^  will  include 
about  95  percent  i~-x  99  percent,  respectively,  of  the  data 
points.  Standard  errors  are  multiplicative  because  of 
the  form  of  the  equation  which  has  been  used: 

PM  =  KR"” 

where 

PM  =  peak  particle  motion  in  g,  cm,  or  cra/sec 
depending  on  typ>e  of  motion 
R  -  distance  in  meters 

K,  n  =  constants  determined  by  regression  analysis 

This  equation  is  linear  in  logarithmic  form 
log  PM  =  log  K  -  n  log  R, 

The  standard  error  of  estimate  is  derived  from  this 
linear  equation  and  is  itself  a  logarithm.  Therefore, 
log  PM  +  log  ay  means  (PM)  (ay)  or  PM/a^. 

The  1  cTy  band  width  is  shown  on  the  graphs  and  the 
value  of  CTy  is  presented  with  the  equations.  In  all 
cases  where  the  regression  equation  is  statistically 
significant,  the  value  of  ay  is  less  than  two.  Thus, 
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tht'  equations  represent,  within  a  factor  of  U’ss  than  two, 
some  68  percent  t,)f  the  data  values  recorded  it  normal  dis¬ 
tribution  is  assuinc'd. 

3.2.7  The  R  a  t  io  of  (  u^ /ad  )  Where  u  j._s  i  c  le 

Ve  1  PC  i  t  y  ,  a  is _ F*ar  tic  1  e  A  cc  e  1  e  ration,  an  ^  d _ ^i.  ^  Par  tic  ^ 

Displacement .  This  ratio  was  examined  for  the  three  sta¬ 
tions  at  which  independent  raeasur oments  of  ihe  three 
ground  motion  parameters  were  made.  These  data  are  listed 
in  Table  3.5.  The  computed  ratios  vary  from  ,076  to  ,68, 
However,  insufficient  data  were  available  to  make  meaning¬ 
ful  statements  about  ratios  associated  with  specific  com~ 
ponents . 

3.2.8  Comparison  of  Observed  Data  with  Predictions, 
Figures  3.16  through  3.18  show  the  comparison  between 
equations  used  for  predicting  Shoal  surface  motions  and 
the  least-square  regression  lines  for  observed  peak  re¬ 
sultant  vectors.  All  prediction  equations  are  higher 
than  the  least-square  fits  for  hardrock.  The  accelera¬ 
tion  and  displacement  predictions  are  nf)  higher  than  the 
+I  Hy  level  and  have  slopes  slightly  in  excess  of  those 
observed.  The  velocity  prediction  is  higher  than  the 
hardr ock-obser ved  peak  resultant  vectors  by  a  factor  of 
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Jn  aJJ  thr(?fc' 


aboMt  two  with  nearly  equivalent  slopes, 
comparisons  the  equations  for  data  observed  on  alluvium 
are  higher  than  the*  prediction  equations  by  factors  in 
the  range  of  1.5  to  4.0. 

On  the  basis  of  this  comparison,  it  is  believed 
that  the  prediction  equations  for  acceleration  and  dis¬ 
placement  adequately  represent  a  reasonable  upper  limit 
of  motion  measured  on  hardrock.  The  prediction  eo’iation 
for  velocity  is  felt  to  be  too  high.  If  the  +1  Cy  level 
is  assumed  to  be  a  reasonable  upper  limit,  then  the  co¬ 
efficient  of  the  velocity  prediction  equation  is  too  high 
by  an  average  factor  of  1.6. 

Based  on  this  data,  the  prediction  equatit^ns  which 
will  be  used  for  future  events  in  granite  will  be  as  fol¬ 
lows  : 

for  points  of  interest  on  hardrock 
a  =  5.03  X  10^ 

d  =  1.51  X  lO^ 

u  =  5.40  X  10<^  W*'^3  r-1.87 

for  points  of  interest  on  alluvium 
a  =  1.06  X  10^  R"2-® 

d  =  5.59  X  10^ 

u  =  1.84  X  10^  R~^*®^ 
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where  units  are: 


a  =  acceleration  in  g 

ci  =  displacement  in  cm 

u  =  particle  velocity  in  cm/st'c 

W  =  yield  in  kt 

R  =  distance  in  meters 

Future  research  and  analysis  may  allow  revision  of  the 
above  equations  so  as  to  include  variables  (e.g.,  frequen¬ 
cy)  whose  effect  is  not  reflected  in  current  prediction 
equations . 
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CHAPTER  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

Detonation  of  the  12.5  kt  Shoal  event  in  granite  has 
provided  much  needed  in/or mati or*  on  seismic  effects  relat¬ 
ing  to  public  safety.  In  addition  to  being  the  best- 
instrumented  shot  to  date  in  granite,  Shoal  also  provided 
the  first  sufficient  amount  of  data  for  comparing  the  am¬ 
plitudes  of  ground  motions  measured  on  hardrock  ’.'ith 
those  measured  on  alluvium. 

The  regression  data  are  considered  to  be  reasonably 
reliable  as  evidenced  by  the  values  obtained  for  the 
standard  error  of  estimates. 

It  appears  that  peak  motions  on  the  alluvium  are  due 
to  the  same  type  of  wave  as  peak  motions  on  hardrock. 

The  equations  used  by  Rolauid  F.  Beers,  Inc.  to  make 
predictions  agree  well  with  observations. 

The  evaluation  of  containment  for  the  Shoal  ev^^nt 
and  the  adequacy  of  the  stemming  plan  were  substantiated 
by  the  results  of  the  test.  Limited  post-shot  exploration 
indicated  that  predictions  of  cavity  radius,  chimney 
height,  depth  of  spalling,  and  radius  of  spalling  were  ac¬ 
curate.  More  extensive  post-shot  exploration  would  permit 
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comparisons  of  other  predicted  effects  which  could  not 
be  evaluated  presently. 

Details  of  data  processing  methods  and  studies  of  a 
long-range  research  nature  will  be  published  in  a  separate 
report. 
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APPENDIX  A 


ANALYSIS  OF  SPALLING  MIt  CHAN  ISM 

Spalling  is  herein  defined  as  the  fracturing  a  solid  by 
tensional  stresses  produced  when  a  seismic  pulse  is  re¬ 
flected  from  a  free  face.  The  conversion  of  a  comf>r  ess  ion - 
al  pulse  to  a  tensional  one  is  considered  to  be  perfectly 
elastic  with  no  loss  of  energy  into  the  atmosphere  at  the 
solid-air  interface. 

The  problem  is  to  analyze  the  stress  situation  at 
continuous  distances  back  from  the  free  face  to  determine 
if  sufficient  tension  exists  to  fracture  the  solid.  The 
tensional  stress  arises  only  from  the  reflected  pulse, 
whereas  those  forces  resisting  fracturing  are  that  portion 
of  the  compress ional  pulse  waveform  which  has  not  yet 
arrived  at  the  interface,  the  tensile  strength  of  the  medi¬ 
um,  and  any  overburden  pressure  which  becomes  significant 
when  horizontal  surfaces  autd  long  wave  lengths  are  con¬ 
sidered  . 

If  the  problem  is  restricted  to  one  involving  com- 
pressional  pulses  arising  from  subsurface  nuclear  detona¬ 
tions  and  a  rock  medium,  generalized  equations  may  be 
used  to  determine  the  stress  condition  at  continuous  depths 
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below  the  surface.  The  equations  are  based  on  hydrodyna¬ 
mic  considerations  and  machine  computations;  in  the  form 
shown  they  are  applicable  at  ranges  of  a  few  hundred  meters. 
The  pulse  shape  assumed  is  one  with  an  exponential  decay 
thus  producing  a  maximum  pressure  the  instant  the  peak 
reaches  a  point  of  interest.  The  stress  condition  at  any 
point  below  the  surface  may  then  be  given  by; 

^net  ~  ^ovb  ^  ^ts  ^  ^pwf  “  ^ref 


where 


^net  -  net  stress  at  the  section  of  interest 

=  constant  -  tensile  strength  of  the  medium 

'^ovb  “  (depth)  -  overburden  stress 
X  2 

^pwf  ~  “  P>^®ssure  wave  form,  or  that  portion 

of  the  pulse  which  has  nc^t  arrived 
at  the  free  suv^ace 

^ref  ~  CR~^'^  -  reflected  peak  pressure 

where 

R  =  depth  of  burial  +  the  depth  of  interest 
r  =  depth  of  burial  -  the  depth  of  interest 
The  value  of  each  of  these  quantities  may  then  be  com¬ 
puted  as  functions  of  depth  to  determine  if  becomes 

>  zero,  in  which  case  a  spall  is  predicted  to  occur.  When 
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a  spall  occurs,  the  overburden  pressure  must  be  reduced 
to  zero  since  the  overlying  strata  are  detached,  and  re¬ 
computed  for  further  depths  of  interest  on  the  basis  of 
the  newly  created  free  surface.  Addi t iona J ] y ,  the  peak 
reflected  pressure  roust  be  reduced  to  the  value  of 

the  pressure  wave  form  at  the  new  free  face  since 

the  energy  in  the  descending  pulse  is  trapped  in  the 
spalled  section.  On  the  basis  of  these  new  values  further 
depths  of  interest  may  be  investigated  to  determine  if  ad¬ 
ditional  spalls  will  occur. 

Since  actual  test  conditions  are  not  similar  to  frce- 
field  predictions  for  a  homogeneous  medium,  corrections 
for  the  variation  in  the  stratigraphic  section  are  intro¬ 
duced.  The  transmission  coefficient,  or  change  of  pres¬ 
sure  across  an  interface,  is  given  by: 

P2  ^°1  ^1 

'^1  ^1^1  *  ^2^2 

where  is  the  density  amd  is  the  compress ional  wave 
propagation  velocity  of  the  original  medium  and  p^,  C2 
are  those  of  the  medium  across  the  interface.  Corrections 
for  transmission  coefficients  must  be  introduced  at  all 
depths  where  they  are  applicable,  to  both  the  ascending 
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and  descending  pulses. 

The  analysis  as  outlined  above  is  confined  to  the  case 
of  normal  incidence  only  and  three  methods  are  presently 
available  for  the  solution  of  the  calculations;  these  are 
by  graphical  procedures,  analog  computer  and  digital  com¬ 
puter  . 


POST-SHOT  CHIMNEY  CALCULATIONS 


Calculations  to  determine  the  height  to  which  post¬ 
shot  chimney  formation  will  occur  are  based  on  the  assump 
tion  that: 

1)  The  post -shot  cavity  is  spherically  symmetrical, 

2)  The  chimney  is  cylindrical  with  a  radius  equal 
to  that  of  the  cavity,  and 

3)  An  increase  in  volume  will  occur  upon  breaking 
of  solid  material. 

The  volume  of  the  cavity,  or  void,  is 


V  =  4/3  tT  R- 
cav  ' 


The  volume  of  the  chimney  is  then  given  by 


''ch  =  ” 


where 


h  =  the  height  above  the  working  point 
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Since  the  chimney  volunoe  includes  one  half  of  the 


cavity  void,  the  volume  of  solids  is  seen  to  be: 


''sol  =  ''ch  -  i  ''cav 

The  original  volume  of  solid  material  must  then  swell 
upon  breaking  to  include  both  its  original  volume  and  the 
volume  of  the  existing  void 


(n)  V 


sol 


cav 


+  V 


sol 


where 


(n)  -  the  swell  factor  for  the  medium  being  con¬ 
sidered  . 

Then  n  j^n  R^h  -  2/3  tt  =  4/3  tt  +  n  R^h  -  2/3  n  R^ 

and  n  (h  -  2/3  R)  =  2/3  R  +  h 

....  h  *  2/3  R  .  .  2/3  ^  (l--n) 

yielding  n  =  iT'T'ITTr  ^  ~  n-1 

The  value  of  (n)  may  then  be  determined  for  mediums 
in  which  experience  is  available  emd  subsequently  be  used 
to  determine  the  chimney  height  for  future  events ,  in  the 
same  medium,  based  on  the  predicted  cavity  radius. 


PERMANENT  DISPLACEMENT  AND  SIRAIN  CALCULATIONS 


With  the  detonation  of  aut  underground  nuclear  explo 
sion  a  post-shot  cavity  results  which  is  larger  than  the 
original  shot  chamber.  The  expansion  of  this  cavity  to 
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its  final  size  entails  the  vaporization,  melting  and 
deformation  of  the  material  surrounding  the  original  cham¬ 
ber.  The  major  portion  of  the  material  is  displaced  radial¬ 
ly  from  the  source.  The  difference  in  position  between  an 
original  point  in  the  medium  and  its  final  location,  as 
measured  upon  re-entry,  is  termed  the  permanent  radial 
displacement  of  the  point. 

EQUATION  DERIVATION 

Given  the  initial  chamber  radius  and  final  cavity 
size,  the  assumption  can  be  made  that  the  medium  has  be¬ 
haved  as  an  incompressible  material  arid  no  loss  of  volume 
of  material  has  occurred  through  vaporization  or  melting. 

The  volume  of  material  which  originally  occupied  the  post- 
shot  cavity  must,  therefore,  be  distributed  throughout  the 
medium.  Assuming  spherical  symmetry  the  following  equation 
may  be  derived: 

Given:  =  radius  of  original  shot  chamber 

R|  =  radius  of  post-shot  cavity 

=  original  radius  to  point  of  interest 

^2  =  final  radius  to  p>oint  of  interest 
If  the  original  and  final  'volumes  are  equal,  then 
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4/3  TT  (R^)^  -  (R^)^  *=  4/3  rr  (R  ^ ^  -  (R|)3 

3  3 

and  (R^)  =  (R2)3  -  (R^)  +  (Rt)3  yielding  AR2  =  ^2  "  ^2 

which  is  the  permanent  radial  displacement  of  a  point 

whose  initial  position  was 

Given  the  predicted  radial  permanent  di sp] acenK?nts  , 

the  tangential  strains  may  be  computed  by  considering  the 

change  in  circumference  at  any  radius  of  interest. 

The  compressive  strain  due  to  a  radial  displacement 

is  obtained  by  plotting  the  displacement  versus  distance, 

and  fitting  an  equation  to  that  segment  of  interest.  At 

distances  relevant  to  radial  fracturing,  this  portion  of 

-1  9 

the  curve  may  be  approximated  by  the  equation  D  =  CR 

dP  _  _2  g 

meters.  The  radial  strain  is  then  "  -1.9  CR  '  meters/ 
meters . 
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APPENDIX  B 


These  data  were  coioputed  for  an  anticii>ated  maximum 
yield  of  14.26  kt  but  are  equally  valid  for  12.5  kt . 

Basic  data; 

Depth  of  burial  -  366  meters 

Cavity  radius  -  2.59  x  10  meters 

4  3 

Cavity  volume  -  7.3  x  10  M 

Fixed  water  content  of  Shoal  granite  -  3%  by  volume 
Melting  point  of  Shoal  grauiite  -  704°  (Reference  18) 
Grams  of  water  per  cubic  meter  of  rock  -  3.0  x  10^ 
Depth  to  top  of  water  table  -  297  meters 
Hydrostatic  pressure  at  W.P.  (366-297)  (.1)  =  6,9  bars 
Density  of  granite  -  po  =  2.67  gm/cm 
The  following  computation  is  made  with  the  assumption 
that  only  the  rock  which  was  melted  would  release  its  fixed 
water . 

Heat  of  fusion  for  granite  -  lOO  cal/gm 
Specific  heat  -  0.2 

Ambient  rock  temperature  assumed  -  20^C 

Sensible  heat  -  136  cal/gm 

Enthalpy  -  9.88  x  10^  ergs/gm  @  700®C 
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Waste  heat  density  -  987.9  x  10^  ergs/gin  (2.67 

gm/ctn^)  =  2,64  x  10^®  ergs/cm^ 

From  Hugoniot  data  for  granite,  this  waste  heat  density 

occurs  at  a  pressure  of  about  340  kilobars.  This  pressure 

has  been  determined  from  Beers,  Inc.  hydrodynamic  computa- 

tions  to  occur  at  a  radius  of  3.42  x  10^  cm. 

Volume  of  rock  melted  -  )  =  167  m^ 

3 

Volume  of  fixed  water  vaporized  s  0.03  (167)  =  5  m^ 

Weight  of  water  vaporized  -  5  x  l.Ox  20^  =  5.0  x  10^  kg 
Specific  volume  of  water  @  704°C  = 

From  steam  tables  for  these  conditions  P  =  0.35  bars 

The  above  pressure  of  0.35  bars  is  less  than  the  hydro¬ 
static  pressure  of  6.9  bars  and  would,  therefore,  be  insuf¬ 
ficient  to  keep  out  water. 

At  some  point,  steam  pressure  behind  the  impervious 
layer  formed  by  the  melted  rock  may  be  sufficient  to  cause 
pieces  of  this  layer  to  breeik  out  into  the  cavity  and  water 
might  enter.  Assuming  that  water  flows  in  at  a  rate  compara¬ 
ble  to  that  in  the  E-1  hole,  the  pressure  would  increase  un¬ 
til  it  reaches  6.9  bars  at  which  time  inflow  would  be  stop|>ed 
by  equalization  of  pressure.  The  time  to  reach  this  pres¬ 
sure  may  be  computed  as  follows : 

At  704°C,  the  specific  volume  at  6.9  bars  absolute  =  0.652  m'^/kg. 


X  10^  ,  .  . 

- ^  =  14.6  m-^/kg 

X  10^ 
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Weight  of  water  =  10~~^  ~  i*12  x  10^  kg 

Additional  water  required  1.12  x  10-*  -  5.0  x  10^  = 

1.07  X  10^  kg 

The  October  4,  1963  PfWO  Construction  Progress  Report 

shows  Hole  E-1  otakiitg  2C3  lit«t/iala  =  265  kg /min, 

1  07  X  10^ 

Time  to  reach  6,9  bars  =s  — * - r  =  404  minutes  m  7  hours 

2.65  X  10^ 
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TABLE  2,1  CX)MPARIS(»J  OF  PHYSICAL  reOPERTIES 
OF  SHOAL  AND  CLIMAX  STOCK  C»ANITES 
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An  average  figure  for  any  typical  granite 


TABLE  2.2  PREDICTIONS  OF  AMPLITUDES  OF  GRajND  MOTION 


Peak  Surface  Particle  Motions  Expected  at  Locations 
of  Interest  ~  Shoal  Event,  12.5  KT 


Location 

Distance 
(Meters ) 

Acceleration 

(9) 

Displacement 

(cm) 

»/elocity 
(cm/sec ) 

Site  Conn.  Rep. 

1,463 

1.38 

3.26 

65.9 

C.P. 

2,240-2,360 

.59-. 53 

1.61-1.47 

29.7-26.9 

Observation  Post 

2,440 

.50 

1.39 

25.3 

Onset  Dauage 

3,900 

-- 

-- 

11.0 

Onset  Daaage 

5,430 

.10 

-  - 

Hwy  50 

7,200-9,600 

.05 7-. 032 

.21-. 13 

3. 4-2.0 

Frenchman's  Sta. 

12,900 

.018 

.078 

1.1 

Nev.  Scheellte 

20,900 

.0068 

.034 

.46 

Salt  Wells 

26,000 

.0044 

.023 

.30 

Fallon  Tower 

35,000 

.0024 

.014 

.17 

Limit  Perceptibility 

54,300 

.OOlO 

-- 

— 
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TABLE  3.1  PEAK  SURFACE  PARTICLE  ACCELERATION  IN  g' 
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TABL2  3.2  PEAK  SURFACE  PARTICLE  D1SPLA(•.EME^3T  IN  CENTIMETERS 
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TABLE  3.4  STATIONS  USED  IN  COMPARISON  OF  DATA  MEASURED  PARALLEL 
AND  PERPENDICULAR  TO  DIRECTION  OF  MAJOR  FAULTING 


Orientation 

Station 

Parallel 

|a 

lA 

1 

2A 

3A 

1 

1 

jC 

1C 

2C 

3~}C 

Perpendicular 

IB 

2B 

1  3B 

4B 

6B 

ID 
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TABLE  3.5  RATIO  OF 

ad 


2 

a  d  u  u 


Station  Comp>onent  Acceleration 

Displacement  Velocity  ad 

2 

ca/sec 

cm  cm/sec 

3B  Z 

77.3 

- 

2.00 

- 

R 

45.2 

- 

2.41 

— 

T 

62.5 

- 

1.79 

— 

vector 

83.3 

— 

3.  12 

4B 

Z 

30.8 

- 

1.27 

- 

R 

36.4 

.232 

2.40 

.68 

T 

30.6 

.124 

1.22 

.39 

vector 

41.3 

.239 

2.47 

.62 

6B  Z 

20.6 

- 

.344 

- 

R 

11.8 

.127 

.418 

.12 

T 

9.72 

.0766 

.354 

.17 

vector 

22.1 

.138 

.480 

.07( 
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Ncvol 
Air  Stotion 


Figure  1.1 

Location  cf  In8ti'u»ent  Station 


Permanent  Rodial  Displacement,  Meters 


Figure  2.1  Predicted  and  Observed  Permanent  Radial 
Dlsplac«nent  and  Predicted  Strain  for  the  Hardhat  Event 
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Radial  Strain,  Percent 


^'•rtnanent  Rodial  Oisplocement,  Meters 


Radiot  Oieltonce,  Meters 

Figure  2.2  Predicted  Strain  and  Permanent 
Radial  Displacement  for  the  Shoal  Fvent 


-  78- 


Radial  Strom ,  Percent 


140 
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Peak  Surface  Particle  Acceleration, 


,  Peak  Surface  Particle  Displacement,  cm 


Figure  2,5  Predicted  Peak  Surface  Particle 
Displacement  as  a  Function  of  Horizontal 
Distance,  Shoal  Event 
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kilometers 

O  25  50 

I  I  I  I  I  I  I 


SCALE 


Figure  2.7  Area  Map  Showing  Predicted 
Distances  to  Significant  Ground  Motions 
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t  P»ok  Surfoce  Porticis  Acceleration, 


10^  10^  10^  10* 


R,  Rodiol  Olstonce,  Meters 

Figure  3.1  Vertical  Component  of  Acceleration 
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Peak  Surface  Porticle  Acceleration, 


R,  Radial  Dittonce,  Meters 

Figure  3.2  Radial  Component  of  Acceleration 


,  Peok  Surfoce  Particle  Accelerot ion, 


R,  Rodiol  Distance,  Meters 

Figure  3.3  Transverse  Component  of  Acceleration 


,  Peak  Surface  Particle  Acceleration, 


10^  10^  10^  10® 

R,  Radial  Oistonce,  Meters 

Figure  3.4  Resultant  Vector  of  Acceleration 
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,  Peak  Surface  Porticle  Displacement, cm 


■o 


R,  Radiol  Distance,  Meters 

Figure  3.5  Vertical  Component  of  Displacement 
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,  Ptok  Surface  Particle  Displacement,  cm 


XI 


R,  Radiol  Distance,  Meters 

Figure  3.6  Radial  Component  of  Displacement 
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P«ok  Suffoce  Porticie  DispSacement 


- -  O  Hardroch 

d  =  I  82*iO^  R"' 

0^=  (.564 

■  ~  ~  A  Al luvium 

d  =  6.32  *  I05  R--'  55 
0^  =  I  621 

±  I  OJj  (Standard  error 
o'  estimate) 


•  -  -  .  • 

I 


I 

I  I  .  I  I 


J . 


R,  Rodiol  Distonce,  Meters 

Figure  3,8  Resultant  Vector  of  Di sr>l acement 


-  •  Hardrock,  Observed  Motion 

O  Hardrock,  Derived  Motion 
li  =  1 . 84  t  10^  R-i  65 

CTy  =1,396 

A  Alluvium,  Derived  Motion 
u  =  I  42  X  10®  R~' 

CTu  =  I.  198 

\ 

^  i:  I  Oil  (Standard  error 
Oi'  estimate ) 


1 1 1 1  ]ili 


10^  )0-» 

R,  Rodioi  Oistonce,  Meters 
Figure  3.9  Vertical  Component  of  Velocity 


Peok  Surface  Particle  Vetocity,  cm/se 


•  Hardrock,  Observed  Motion 

OHordrock,  Derived  Motion 
u  =  3.38  X  10®  R“' 

Cu  =  1.478 

^  Alluvium,  Derived  Motion 
u  =  6.21  X  10^  R"' " 

CTy  =  I  143 

^  (Stondord  error 
of  estimote) 


imr 

I  III 


-H-l-n - 

l+ll-b- 


r  f  4  -t 

I  I  >  I 


1  .IJ  LiiJ 


R,  Rodiol  Distance,  Meters 

3.11  Tiansversp  Component  of 
Veloci ty 


Peok  Surfoce  Particle  Velocity,  cm/s 


,  Peok  Surfoce  Porticle  Acceleration, 


R,  Rodiol  Dislonce,  Meters 

Figure  3.13  Comparison  of  Resultant  Vector  of  Acc«’ 1 er a t ion 
Measured  Parallel  and  Perpendicular  to  Major  Fault  Trend 
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,  Peok  Surfoce  Particle  Displacement,  c 
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,  Peok  Surface  Porticie  Velocity,  cm/sec 


Figure  3.15  Comparison  of  Resultant  Vector  of  V'elocity 
Measured  Parallel  and  Perpendicular 
to  Major  Fault  Trend 
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,  P«Qk  Surface  Particle  Accelerotian, 


Figure  3.16  Comparison  of  Predicted  and  Observed 
Resultant  Vector  of  Acceleration 
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,  Peak  Surface  Particle  Displacement, 


Figure  3.17  Comparison  of  Predicted  and  Observed 
Resultant  Vector  of  Displacement 


lOO 


,  Peak  Surface 
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